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1. INTRODUCTION 

Due to the individual merits such as no contact, no wear, no lubrication and adjustable dynamics, 
magnetic bearings (MBs) are being widely applied to a variety of industrial products, especially the high- 
speed rotating machinery and precise positioning systems [1-3]. Considering the inherent nonlinearity and 
open-loop instability of MBs, the accurate modeling and control are two essential aspects of achieving high- 
performance requirements. Nowadays, PID controllers have received considerable attention in the last years 
both from an academic and industrial point of view [1-5]. In fact, in principle, they provide more flexibility 
in the controller design, concerning the standard PID controllers, because they have five parameters to select. 
However, this also implies that the tuning of the controller can be much more complicated. They have been 
successfully applied in practical applications such as motion control of manipulators and chaos control of 
electrical circuits. In these applications, it has been verified that PID controllers can improve the performance 
of traditional control system adopting integer order PID controllers. The most important advantage of the PID 
controllers is that they can afford broader possibilities offered by their new fractional order dynamics [4-7]. 
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However, this also indicates that the tuning strategies of PID controllers are much more complicated. In the 
researches on the PID controllers, tuning of controller parameters has become a significant issue. 

In general, the tuning methods for PID controllers are classified into analytical, numerical, and rule- 
based ones. In [6-7] the controller parameters have been analytically derived by solving nonlinear equations 
fulfilling the gain/phase crossover frequency and phase/gain margin specifications. The robustness to loop 
gain variations specification proposed in [8] has also been widely used to design PID and proportional— 
integral (PI) controllers. The merits of the analytical method are apparent; however, it is available only when 
the equations are few in number and simple. Therefore, it is challenging to obtain a complete PID controller 
for the MB system by solving five complicated nonlinear equations. As for the rule-based method, it can 
easily calculate the controller parameters based on empirical tuning rules, which can be observed in [9-21]. In 
this paper, we propose and investigate the optimal tuning PID controller of unstable fractional order system 
by desired transient characteristics using the real interpolation method (RIM). The main advantages of this 
method are drawn as the followings: 

1) Carrying out an investigation of the stable region of coefficients of a PID controller using D- 
decomposition method. 

2) Applying the method to investigate an unstable fractional order system. 

The rest of this paper is organized as follows. In Section 2 the problem is formulated. Numerical 
results and discussions are described in Section 3. Conclusions are drawn in Section 4. 


2. PROBLEM FORMULATION 
In this section, Figure | illustrates the control system with the negative unity feedback. The 
mathematical description of this model is presented. 
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Figure 1. Control system with negative unity feedback 


The transfer function of the PID controller can be formulated as 
K, 

C(s)=K,+—+K,s (1) 
sS 


The transfer function of the plant (in the fractional form) is calculated as 


_ BGG) © bs +b_s +---+b5% +b,s” 
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m m-1 


G(s) 





(2) 





where m,n E N,and ap, ..., am bo, --, Dn: m > am-1 > > Ay > Qo = 0, Pn > Pn-1 > 0 > Ba > Bo = 
0 are arbitrary real numbers am # 0, b, # 0. 
The transfer function of PID controller (1) can be rewritten in the rational form as the following 


B(s)_ K,s°+K,s+K, 





C(s) = = 3 
Ora ; (3) 
The characteristic equation in this analysis can be expressed as 

W(s) =1+ C(s) G(s) (4) 


Characteristic polynomial is formulated as 





Determination of all stability region of unstable fractional order system in case of ... (Van-Duc Phan) 


1736 O ISSN: 2088-8694 


P(s) =s4,(s)+(K, s°+K, s+K,)B,(s) (5) 


A gain-phase margin tester (GPMT) can be thought of as a “virtual compensator”, provides 
information for plotting the boundaries of constant gain margin and phase margin in a parameter plane [14]. 
The frequency independent GPMT is given as in [15]: 


GM, $) = Me* . (6) 


For a given IOPID controller parameters K ,,K;, K, the closed-loop system is said to be bounded- 
input bounded-output (BIBO) stable if the quasi-polynomial P(s,K ,,K;,K,,) has no roots in the closed right- 
half of the s-plane (RHP). The stability domain S in the parameter space P with K,,K,,K, being 
K,,K,) all lie in 


open left-half of the s-plane (LHP). The boundaries of the stability domain S which are described by real root 
boundary (RRB), infinite root boundary (IRB) and complex root boundary (CRB) can be determined by the 
D-decomposition method [15], [16] . These boundaries are defined by the equations P(0,K) =0, P(#,K) =0 


and P(+ja@,K)=0 for@e(0,%), respectively, where P(s,K) is the characteristic function of the closed- 


coordinates is the region that for K,,K,,K, €S the roots of quasi-polynomial P(s,K,, 


loop system and K the vector of controller parameters. 


2.1. Determining RRB 
In applying the descriptions of stability boundaries of the stability domain S to the FOCE in (1), the 
RRB turns out to be simply a straight line given by 


P(0,K ,K,,K,)=0 = K, =0, (7) 


p? 
for s“ =1 in the transfer function of the plant in (2). 


2.2. Determining IRB 

The calculating of the IRB faces with more technical difficulties due to the fractional component. 
FOCE possesses an infinite number of roots, which cannot be calculated analytically in the general case. 
However, the asymptotic location of roots far from the origin is well known [17], which may lead to IRB. 
The objective of this section is to determine the stabilizing region in (Kp, K;) plane with given Kg, and values 
for which the following complex polynomial is Hurwitz: 


D(s) = sL(s)4,(s)+(K, s*+K, s+K,)M(s)B,(s) (8) 


where L(s) and M(s) are given complex fractional order polynomials. When L(s) = M(s) = 1, the 
stabilization of Eq. (8) reduces to the standard IOPID stabilization. 


From (8), we have D(%) =0. Suppose whens > œ, [ L(s) 4, (s) | / [ M(s) Be (s) | — cs‘ where t 
and c are real and complex numbers respectively. We have [16]: 

(a) If t #1, then the boundary does not exist. 

(b) If ¢=1 and c is not real, then the boundary does not exist. 

(c) If ¢=1 and c is real, then 





N(s)(K,s+K,+K,s° 
lim Dis) =liml4 ORAK EK E (9) 
s>% sL (s) A, (s) s25 SAG (s) 
kas (10) 
c 


The RRB and IRB lines can be obtained from (7) and (10), respectively. 


2.3. Determining CRB 
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To construct the CRB, we substitute s = j@ into (5) to obtain 


P§a) = 4.9 @)4,G.@) + B-G@)B, Ga) =0 (11) 

Using D-decomposition to find stability [15]. In the s plane j’ "i = cos(5-9) + isin 9) ; 
where @is a real number. The formula (6) can be expressed into a formula with the separate real and 
imaginary components 

Pia) =| (K,-K, ©) Ryg(@) - @1,¢(@) |+ fo] K,Ryg(@)+Ryg(@)|=0, A2 
where R,g  ,Ryg,lyg are the real component of B,(j@), the real and imaginary components of 
A, (j@) , respectively. 


Finally, by setting the real and imaginary parts equalized zero the formula (12) leads to the 
detail form: 


te (@)+K,R,_(@) =0 
(13) 


(K,-K, @ )R yg (@) -@1,,(@) =0 


The solution in the (K;, Ki) plane 


This is a three-dimensional system in terms of the controller parameters K |, K,, K, . In order to deal 


p? 


with two unknowns K.,, K,, the value of K, will be fixed to find the stability region in the (K „, K,) plane. 
Equation (11) rewrite in terms of two unknowns K „, K, . 
K, =-R 4¢(@)/R5¢(@) (14) 
K, = @1,¢(@)/Ryg(@)+K, @° 


The above two equations trace out a curve in the (K_,,K,) -plane representing the CRB, for fixed 


p? 
K,, as @ run from 0 too. 

The solution in the (Kp, Ka) plane 

Next, to find the stability region in the (K;, Ka) plane, the value of K; will be fixed. Then formula 
(13) can be rewritten as: 


: = -R yg(@)/ Ryo) (15) 


K, =[K,-@1,¢(@)/ Ry (@) |/ 


The above two equations trace out a curve in the (K,, K,) -plane representing the CRB, for fixed 


K,;, as @ run from 0 too. 


The solution in the (Ki, Ka) plane 
Next, to find the stability region in the (Kj, Kz) plane, the value of Ka will be fixed. Then formula 
(13) can be rewritten as: 


i = -R (0) /Rpo(0) (16) 


K, =[K,-@1j¢(@)/Ryg(@) |/ 


The above two equations trace out a curve in the (K,, K,,) -plane representing the CRB, for fixed 
K,,as @run from 0 tom. 


ar 


3. NUMERICAL SAMPLE AND DISCUSSION 
Given fractional order transfer function of an unstable bearing system [18] [19], 
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6438. 
G(s) = 17 
6s z +330.04s'*! —84268.855°" —15869154.14 m 
Investigation of the stability region of the unstable fractional order transfer function of the bearing system 


(17) is carried out in the following parts. 





3.1. Determining RRB 
It is important to determine the stability region of the PID controller, tuned for the system (17). The 
RRB can be determined by (7) as the bellow 


K,=0. (18) 


3.2. Determining IRB 
The IRB can be determined according to (9) 





643.8(K s+K. +K s 
lim ita (K,s+K,+K,s") 


=1 19 
s+» sD(s) s> (577° + 330.04s'*! —84268.855°" ~15869154.14)) (19) 


Not exist the IRB. 


Ki 
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Figure 2. The CRB graph in the (Kp, Kj) plane and Kg = 10. 


3.3. Determining CRB 

The CRB can be determined by three conditions (13), (14) and (15). Figure 2 shows the stability 
region of the transfer function (17). The CRB in the (K_p,K_i) plane when K_d=10 is shown in the Figure 2. 
The cross-line area is a feasible region for the coefficients of the PID controller. The cross-line region 
demonstrates stability region of K_p,K_i where K_d=10. Obviously, the stability region of PID controller’s 
coefficients lie in the first quarter of the coordinate [22-25]. 

Furthermore, Figure 3 shows a more feasible area with five values of K4. 




















Int J Pow Elec & Dri Syst Vol. 10, No. 4, Dec 2019: 1734-1741 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1739 





Figure 3. The CRB graphs in the (Kp, K;) plane with 5 values of coefficients Ka. 
Figure 3 shows that by increasing the value of the coefficient Kq , the feasible region is expanded. 





Figure 4. The CRB of the PID coefficients 


In Figure 4. The distribution of K,, K;, Kq in the 3-D system is represented. In more details, we can 
investigate the stability region of the PID controller’s coefficient separately in the 2-D coordinate plane. In 
the Figure 4. 





Figure 5. The CRB in the ce Ki) plane 


Figure 5 illustrates the stability region of K, and K;. Obviously, the stability region is covered the first 
and fourth quarters of the coordinate system. From the Figure 5, there is a significance that the coefficient 
K, >= 2465. Furthermore, The Figure 6 show the stability region in the Kp, Ka plane. The values of the 
coefficient K,,, Kq are located in the first and fourth quarters, and of course K, 2 2456. 


50 





Figure 6. The CRB in the (Kp Ka) plane 





Determination of all stability region of unstable fractional order system in case of ... (Van-Duc Phan) 


1740 


g ISSN: 2088-8694 


The CRB in the (K;, Kg) plane shows that the distribution of K; and Kg 


4. 





Ki 


Figure 7. The CRB in the (K;, Ka) plane 
The results simply lead to the conclusion 


K, >= 2465 
| ” (20) 


K,>0 


Three parameters K,,K,,K, of the PID controller are determined in four figure 4-7 [22-25]. 


CONCLUSION 


In this paper, we propose and investigate optimal tuning PID controller of unstable fractional order 


system by desired transient characteristics using the real interpolation method (RIM). The research shows 
that the main advantages of this method are drawn as the followings: 1) Carrying out an investigation of the 
stable region of coefficients of a PID controller using D-decomposition method; 2) Applying the method to 
investigate an unstable fractional order system. 
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